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ABSTRACT 

Aims. We present the performance of the Integral Field Spectrograph (IFS) of SPHERE, the high-contrast imager for the ESO VLT 
telescope designed to perform imaging and spectroscopy of extrasolar planets, obtained from tests performed at the Institute de 
Planetologie et d’Astrophysique de Grenoble facility during the integration phase of the instrument. 

Methods. The tests were performed using the instrument software purposely prepared for SPHERE. The output data were reduced 
applying the SPHERE data reduction and handling software, adding an improved spectral deconvolution procedure. To this aim, we 
prepared an alternative procedure for the spectral subtraction exploiting the principal components analysis algorithm. Moreover, a 
simulated angular differential imaging procedure was also implemented to estimate how the instrument performed once this proce¬ 
dure was applied at telescope. The capability of the IFS to faithfully retrieve the spectra of the detected faint companions was also 
considered. 

Results. We found that the application of the updated version of the spectral deconvolution procedure alone, when the algorithm 
throughput is properly taken into account, gives us a 5cr limiting contrast of the order of 5x10“'’ or slightly better. The further appli¬ 
cation of the angular differential imaging procedure on these data should allow us to improve the contrast by one order of magnitude 
down to around 7x10“^ at a separation of 0.3 arcsec. The application of a principal components analysis procedure that simultane¬ 
ously uses spectral and angular data gives comparable results. Finally, we found that the reproducibility of the spectra of the detected 
faint companions is greatly improved when angular differential imaging is applied in addition to the spectral deconvolution. 

Key words. Instrumentation:spectrographs - Methods: data analysis - Techniques: imaging spectroscopy - Stars: planetary systems 


1. Introduction 

High-contrast direct imaging of extrasolar planets is very chal¬ 
lenging because of the very high brightness contrast between the 
planets and their host (10 ® for young giant planets with intrinsic 
emission observed in the near infrared and 10“^ -10“^ for old gi¬ 
ant and rocky planets seen through stellar reflected light), and the 
very small angular separation between the planets and their hosts 
(few tenths of arcsec for planet at ~10 AU from their hosts and 
tens of parsec away). For this reason, only 47 low-mass compan¬ 
ions out of more than 1000 planets discovered so far have been 
detected by exploiting high-contrast imaging (see Extrasolar 


Enciclopaedi^B - ISchneider et al.ll201 ll) . A new generation of 
instruments devoted to exploiting this technique is now becom¬ 
ing operational. They combine extreme adaptive optics (XAO), 
providing high Strehl ratio (SR) by correcting high-order wave- 
front aberrations, with high-efficiency coronagraphs to attenuate 
the diffracted light from the host star. To fully exploit the poten¬ 
tial of these instruments, differential ima ging techniques, suc h 
as angular differential imaging (ADI; see iMarois et al.ll2006bl) . 
simultaneous spec tral differential imagi ng (S-SDI; see, e.g., 
iRacine et al .119991 and iMarois et al.l2005h and spectral deconvo- 

' http://exoplanet.eu 
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lutioiJ3(SD; see lSoarks & Fordll2002l and lThatte e t al.||200 7[) are 
also u sed. One of these instruments. Project 1640 dCrenn et al.l 
IMIh . is already operating at the 5 m Palomar Telescope an d 
led to good scientific results (see, e.g.-IOpDenheimer et ^1.1201^. 
The Gemini Planet Imager (GPI; see iMacintosh et alJ 2006h is 
starting science operations at the Gemini South Telescope while 
SPHERE is currently in its science verification phase at Paranal. 
A fourth instrument called CHARJS is expected to beco me op¬ 
erative by the end of 2015 (iPeters-Limbach et al.ll2013h at the 
Subaru Telescope. This instrument is expected to expand the ca¬ 
pabilities o f SCExAO, alread y available at the same telescope 
(see e.g. - ICurrie et~aD l2014t) . Besides these instruments, the 
LBT telescope with its lower order adaptive optics system has 
proven to be a very po werful instrument for the imaging of extr a- 
solar planets (see e.g. lSkemer et al.ll2 012l:lEsposito et al]l2013h . 

The instrument SPHERE (iBeuzit et al. |2006|) , dedicated 
to high-contrast imaging of faint planetary mass companions 
around nearby young stars, is installed on the Nasmyth platform 
of the ESO Very Large Telescope (VET). It is divided into four 
subsystems; 

- the Common Path and Infrastructure (CPI) receives and 
transfers telescope light, feeds the other three subsystems 
with a highly stabilized adaptive optics (AO) corrected and 
coronagraphic beam. The Extreme AO system for SPHERE, 
calle d SAXO, produce s a beam with a Strehl ratio better than 
90% (iPetit et ail. 20121). T he AO deformable mirror (DM) is 
described in iHugot et al.l (l2012h . 

- the Zurich IMaging POLarimeter (ZIMPOL) produces direct 
and differential polarimetric images in the visible between 
500 and 900 nm (iThalmann et al.i2008h . 

- the InfraRed Dual-band Imager and Spectrograph (IRDIS) 
has various observing mo des: classical imag ing (Cl), dual¬ 
band imaging (D BI; see IVigan et al.l 20T^ , polarimetric 
imaging (DPI; see |LMiglois_etal[ 20K), or long-slit spec¬ 
troscopy (LSS - see lVigan et al.ll2008l) . It works in the near 
infrared between 0.95 a nd 2.32 /rm with a field of view 
(EOV) of 12.5x11 arcsec dPohlen et al.ll2008ll . 

- the IPS also operates in the near infrared between 0.95 and 
1.35 //m with a 2-pixel spectral resolution of 50 for the YJ 
mode and between 0.95 and 1.65 urn with a spectral res olu- 
tion of around 30 for the YH mode (IClaudi et al.ll200^ . Its 
EOV is of 1.7 arcsec square. The integral field unit (IPU) is 
composed of a double system of lenslets set on an a-focal 
configuration and tailored to reduce cro ss-talk between bot h 
adjacent lenslets and adjacent spectra dAntichi et al.ll200^ . 
The raw data for IPS is composed of around 21000 spectra 
aligned with the detector columns over a hexagonal grid ro¬ 
tated by 10.9° with respect to the disper sion axis. The data 
reduction and handling (DRH) software dPavlov et al.ll2008l) 
transforms each detector frame into a datacube composed of 
39 monochromatic images of 291x291 pixels for both the 
observing modes. 


Among the other modes of SPHERE, the IFS achieves very high 
contrast by combining spectral and angular differential imaging. 
In recent years several teams have dedicated efforts toward de¬ 
veloping spectral extracting tools for high-contrast imaging with 
integral field spectrographs data. For example, ISeifahrt et al.l 

^ The authors are aware that the definition of spectral deconvolution 
is probably inadequate since no deconvolution is actually performed. 
However, this name is routinely used in the high-contrast imaging com¬ 
munity. Moreover, the method is clearly explained in Section [TT] For 
these reasons, we have decided to keep it in our paper. 


(I2nn7h and lLavigne et al.l (l2009h presented an anal ysis aimed at 
extrac ting the spectrum of GQ Lup b. Moreover. iBowler et al.l 
(120101) and iBarman et alJ (1201 Ih performed high-contrast spec¬ 
troscopy of HR8799b. More recently, the spectroscopy of all 
HR87 99 system planets was performed by iQppenheimer et aH 
(12013h and a library of spectr a of 15 young M6-L0 dwarfs was 
compiled by iBonnefoy et al.l (l2014l) . A method for the spec¬ 
tral characterization of faint co mpanions exp l oiting the IFS of 
Projects 1640 was described in Pueyo et ^ (120141) while the 
TLOCI algorithm (iMarois et al.ll2014 ) was deyeloped to exploit 
GPI IFS data, using both the ADI and the S-SDI. A method for 
high-contrast imaging with AD I technique was dey eloped in the 
context of the SEEDS pipeline dBrandt et al.ll2013]) . 

In this paper we present the results of the laboratory tests per¬ 
formed and highlight the performance limits of the IFS in terms 
of contrast and signal to noise (S/N) as determined from lab¬ 
oratory data. We obtained these results with dedicated speckle 
subtraction algorithms, which we will describe. While on-sky 
data from SPHERE are currently becoming ayailable we think 
that this paper exploiting the laboratory data is still useful giyen 
that it allows us to present the state-of-the-art data reduction al¬ 
gorithms for SPHERE IFS. 

The outline of this paper is: in Sectionl^we describe the test 
procedures, in Section[2we describe the data reduction pipeline, 
in Section |4] and Section |5] we present the results of the data 
reduction for the YJ- and YH-mode, and in Section |6] we giye 
our conclusions. 


2. The IFS tests description 

The whole SPHERE instrument is at present in its science yerifi- 
cation phase at Paranal. Before shipment to ESO, it was assem¬ 
bled in a large dedicated clean room at Institute de Planetologie 
et d’Astrophysiquede Grenoble (IPAG). A large numberof func¬ 
tional and performance tests were performed on the instrument 
from the end of 2011 to the end of 2013. The aim of this paper is 
to describe the tests performed on IFS. The first goal of these 
tests was to check that the instrument control software (INS) 
worked properly. It is organized in different components, defined 
as ’’templates”, which can be run autonomously or within ob- 
serying blocks, according to the ope rational model adopted at 
VET (see e.g. iBaruffolo et al.l 1201^ . These templates execute 
all the tasks needed to perform the calibration of the instrument 
and acquire scientific obseryations. 

The main calibration templates perform the following tasks: 
detector dark and flat field acquisition, determination of the spec¬ 
tra positions, wayelength calibration and finally acquisition of 
the IFU flat field. The IFU flat field is required to calibrate dif¬ 
ferent outputs from different lenslets composing the lenslet ar¬ 
ray. Other calibration templates were executed to monitor the 
instrument such as the gain and readout noise (RON), ghosts, 
distortion and background. All these templates were intensiyely 
and repeatedly tested and, when necessary, modifications were 
made to improye their outcome or to correct errors during their 
execution. Other calibration templates dedicated to the obserya¬ 
tions (e.g., sky flat) will be fully tested when the instrument is 
installed at the telescope. 

Two templates enable the acquisition of scientific data: one 
of them is for obseryations of single objects with IFS in the YJ 
mode IRDIFS, while the second is for the obseryations in the 
YH mode (IRDIFS_EXT). All these templates run IFS in parallel 
with IRDIS. Seyeral dichroic beamsplitters in the CPI, located 
after the coronagraphic focal mask and before the Lyot stops, al¬ 
low for the feeding of both instruments simultaneously. IRDIS is 
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Fig. 1. Example of IFS coronagraphic raw data. The displayed 
image covers a square 1.77x1.77 arcsec FOV. This observation 
was performed in YJ-mode. 



Fig. 2. Zoom on the central part of the image from Figure [1] 
showing the spectra more clearly. 


used both for acquisition purposes, like hne coronagraphic cen¬ 
tering and the noncommon path aberrations calibration, as well 
as for scientihc purposes. 

The laboratory tests performed at IPAG used a telescope sim¬ 
ulator (TSIM) including rotating phase screens to emulate atmo¬ 
spheric turbulence corresponding to 0.6", 0.8", and 1.2"seeing 
conditions. The phase screens have adjustable rotation speed to 
simulate wind speeds ranging from 2 to 10 m/s. Stars of differ¬ 
ent brightness were emulated by inserting different neutral den¬ 
sity biters. It is important to say that this method just simulates 
a fainter star from the IFS point of view, while this is not true 
for the AO. This is important given that the AO performance de¬ 
grades for faint stars. Moreover, several coronagraphs (SPHERE 
allows the use of both apodized Lyot and 4-quadrant corona- 
graphs; see iBoccaletti et aP l2008h have been tested. However, 
we report on results obtained only with the apodized Lyot coro- 
nagraph with a mask diameter of 185 mas (projected on sky) 
since the performance of the 4-quadrant coronagraph is not yet 
optimized. Possible explanations for performances loss, com¬ 
pared to expectations, with this type of coronagraph could be 
related to the presence of low-order aberrations, probably a de¬ 
focus, or a possible aging of the optical components materials. 
These reasons could explain the higher level of performances ob- 
tained during the prel iminary tests performed at LESIA (see e.g. 
iBoccaletti et al.ll2008l) . Figure [1] shows an example of IFS (in the 
YJ-mode) coronagraphic raw data, while Figure|2]shows a zoom 
on the central part of the same image highlighting spectral lay¬ 
out. In this case, we used a simulated star with a J magnitude 
of 2.6 with typical seeing and wind speed condition for Paranal 
(0.8 arcsec seeing and 5 m/s wind speed). 

3. Data reduction pipeline 

We performed the data reduction using the DRH pipeline devel¬ 
oped by the SPHERE consortium. This software is composed by 
different modules, called recipes, aimed to reduce specibc data 


sets generated by the SPHERE instrument software. The data 
reduction pipeline includes the following steps: 

- creation of the master dark; 

- creation of the detector master bat beld; 

- debnition of the IFS spectra positions. At this step, each de¬ 
tector pixel is attached to a particular spectrum and a brst 
guess wavelength is assigned to it, based on a model. The 
corresponding data set is obtained by illuminating the IFS 
with a broadband light source connected to an integrating 
sphere located in the calibration arm; 

- wavelength calibration of each spectra. At this stage the 
wavelength assigned to each pixel is rebned using data ob¬ 
tained illuminating the IFS with several laser light sources 
(three for the YJ-mode and four for the YH-mode). Those 
lasers feed the same integrating sphere. 

- creation of IFU bat beld to take the response to a uniform 
illumination of all the IFU lenslets properly into account; 

- reduction of the scientihc images with the creation of 
monochromatic images datacubes. The datacubes are ob¬ 
tained by interpolating the individual spectra spectrally and 
spatially. A more detailed description of the procedure used 
in the DRH to create the calibrated datacube can be found 
in lP^lov et al.l (l2008h : 

- subtraction of the speckle pattern by exploiting the SD 
method. 

To improve the operation of the instrument and the accuracy 
of the calibrations, we requested upgrades of the DRH recipes 
in recent months. Starting from the basis offered by the DRH 
pipeline, we also improved the spectral deconvolution procedure 
with the aim to replicate th e results of the IFS performance sim¬ 
ulations dMesa et al.l201 ih and to improve the overall instrument 
performances. In particular, we improved the algorithms for the 
debnition of the star center, the wavelength rescaling factor, and 
of the full width at the half maximum (FWHM) of the speckles at 
different wavelengths. In addition, the principal component anal¬ 
ysis (PCA) technique has been adapted for the IFS. Both these 
procedures were developed using the IDL language. Their inclu- 
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sion in the DRH pipeline is foreseen in the next years. In the 
following we describe these procedures in detail. We would like 
to emphasize that the procedures described here are the results 
of a series of different attempts aimed to optimize their perfor¬ 
mance. In the following we will describe the current versions of 
the procedures. 

3.1. Spectral deconvolution for IFS 

The aim of the SD algorithm is to reduce the speckle n oise. 
This method is des cribed in detail in I Sparks & FordI (l2002h and 
iThatte et~an(l2007h . It can be summarized by the following steps: 

- each single frame in the datacube is rescaled according to its 
wavelength in such a way that the speckle pattern for the star 
becomes similar at all wavelenghts while its companion is 
located at different positions for different frames; 

- a polynomial fit is performed along the wavelength direction 
for each pixel of the rescaled datacube; 

- the obtained fit is then subtracted from the rescaled datacube 
to reduce the speckle noise; and each frame of the datacube 
is rescaled back to the original dimensions in such a way that 
the companion image is set back to the same position in all 
the frames. 

The main difficulties of the SD method for IFS is to properly es¬ 
timate the center of each monochromatic image in the datacube 
and to rescale the images because the rescaling factor slightly 
varies with wavelength due to imperfect wavelength calibration. 
To overcome these problems, satellite spots were artificially in¬ 
troduced by generatin g periodic phase offset s on the deformable 
mirror as proposed in Ikanglois et alJ (l2013h for what concerns 
their use for SPHERE obse rvations. The use of satellite spots 
was pr eviously proposed by ISi varam akrishnan & Qppenheime^ 
(|200^ and by iMarois et a^l.r(f2006ah . These satellite images, 
symmetric with respect to the central star as shown in Figured 
are used to compute the star center position very accurately (with 
a precision of some tenths of spatial pixel) and, moreover, to 
calibrate the rescaling factor more accurately by exploiting their 
shifting outward with increasing wavelengths. 

A second order of problems consists in the fact that the im¬ 
age could not be uniformely illuminated because of an out of 
focus ghost on the telescope simulator optics. To solve this prob¬ 
lem, we performed a bidimensional fit on the original image. The 
image resulting from this procedure is an inclined plane func¬ 
tion that is subtracted to the original image. While this problem 
arises from the telescope simulator, we cannot exclude that a 
similar problem could arise for the on-sky data. For this reason 
this procedure should be maintained in the pipeline and used just 
in case it is needed. Besides, when rescaled properly according 
to their wavelength, the speckle FWHM varies with wavelength 
because it is not uniquely determined by diffraction. This results 
in an inaccurate subtraction of the speckles. To solve this prob¬ 
lem, all the images are smoothed so that the speckle FWHM re¬ 
mains unchanged for all wavelengths. This method enables more 
precise speckle subtraction, but the smoothing slightly degrades 
the spatial resolution to about 1.5 A/D. Finally, to compensate 
the flux variations within original images and not properly cor¬ 
rected by the subtraction of the inclined plane function described 
above, each image is divided in eight radial sectors for which 
normalization factors are calculated using the median of the im¬ 
age for separation ranging from 6 to 15 A/D. The images were 
then divided by these normalization factors taking care to main¬ 
tain comparable flux for all the images at different wavelengths 
before applying the SD procedure. These procedures generate 



-579 -132 318 765 1215 1662 2108 2559 3005 

Counts 


Fig. 3. Example of a monochromatic (A = 1.047/im) reduced 
image where the four artificial satellite spots are clearly visible. 
This image was taken exploiting the SPHERE apodized Lyot 
coronagraph (ALC). The diameter of the coronagraphic mask 
used is 185 mas. A bar indicating the length of 0.1 arcsec has 
been superimposed to the image to make the scale clearer. 


images more uniform in flux and better aligned with the speckle 
pattern. As a final step, we developed the SD procedure to ex¬ 
clude the planet signal when fitting the function to be subtracted 
from spectra to reduce the speckle noise. This was done by ex¬ 
cluding all the pixels around the maximum of the spectrum from 
the spectrum to be fitted if larger than its median plus fifteen 
times the standard deviation of the remaining part of the spec¬ 
trum itself. The number of excluded pixels is calculated accord¬ 
ing to the fraction of the spectrum that would be covered by a 
possible planet at the separatio n under cons i derati on according 
to the theory of the SD given in iThatte et al.l (l2007h . 

3.2. The PCA for IFS 

As an alternative to the SD algorithm described above, we im¬ 
plemented a P CA procedure to be app lied to the IPS data as de¬ 
scribed, e.g., in ISoummer et al.l (l2012ll . In this case, the data pre¬ 
reduction is identical to what has been described in Section im 
but the normalization of the images is performed using the im¬ 
age median between 6 and 15 A/D, without division in sectors. 
A stellar profile, computed from the median of all the pixels 
equidistant from the central star, is then subtracted from each 
image. Only the central part of these resulting images is consid¬ 
ered to All the two-dimensional array needed to feed the PCA 
procedure. The size of these subimages can be selected by the 
user, but from our experience we have found that the optimal 
size is 200 pixels centered on the center of the star. This allows 
us to work on the most important part of the image near the cen¬ 
ter. Moreover, discarding the most external part of the image, 
we are able to avoid a procedure that is too expensive computa¬ 
tionally. The PCA, relying on the singular value decomposition 
algorithm, generates three arrays which, combined in the correct 
way, lead to the eigenvectors and the eigenvalues used to recon¬ 
struct the original data. The principal components for this recon¬ 
struction are given by the matrix product of these two arrays. A 
principal components subset is used to generate an image with 
the quasi static noise pattern that can then be subtracted from 
the original image. The number of principal components used 
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is adjustable by the user. These adjustments will be discussed 
in more detail in Section 14.21 Clearly, the larger the number of 
principal components used, the better the noise subtraction, but, 
on the other hand, the larger the cancellation of the signal from 
possible faint companion objects. 

3.3. Angular differential imaging 

During the laboratory tests continuous FOV rotation by the dero¬ 
tator was not possible because this would result in a misalign¬ 
ment of the Lyot stop when using the TSIM since the telescope 
pupil is inside the TSIM itself. So, strictly speaking, no imple¬ 
mentation of the ADI technique was possible. We tried to evalu¬ 
ate the efficiency of this algorithm, however, by subtracting im¬ 
ages taken at different epochs. To this aim, we implemented a 
very simple ADI procedure that subtracts from each image all 
the other images in the datacube excluding just the adjacent im¬ 
ages. Each image was normalized before the subtraction divid¬ 
ing for the value of the total flux of the image itself. The analysis 
was performed using the data reduced with the speckle subtrac¬ 
tion methods described in Section 13.11 and 13.21 In this config¬ 
uration, SPHERE was proven to be extremely stable, with no 
significant decorrelation of the instrumental speckle pattern over 
several tens of minutes. Given the great stability of the speckle 
pattern in the laboratory images, which will not be reproduced 
for the on-sky data, these results give an optimistic estimation 
of the on-sky instrument performances, where we expect more 
rapid decorrelation of instrumental speckle s with timescale of 
about ten minutes (see iMartinez et ai]l2013h with derotation. It 
is important to stress that the temporal stability will be very im¬ 
portant for on-sky operation. Given that SPHERE was carefully 
designed with this aim in mind, however, we expect very good 
stability even for on-sky data. 

3.4. Algorithm throughput 

To properly evaluate the contrast level achievable by exploiting 
the speckle subtraction method described above we have to ac¬ 
count for the algorithm throughput. This effect is inherent to the 
differential techniques used in high-contrast imaging. In fact, 
these techniques require the subtraction of two images: in one 
image the faint companion signal is present, and in the other it is 
either absent or weak. The two images should have a noise pat¬ 
tern, mainly due to speckles, as similar as possible to each other. 
This leads to select images with small differences in time for an¬ 
gular differential imaging or in wavelength for spectral differen¬ 
tial imaging. The closer the images are in time and wavelength, 
the more likely the faint companion signal is present in the com¬ 
parison image. Hence, once the two images are subtracted from 
each other, part of the signal from the faint companion is can¬ 
celed together with noise. To optimize the detection, the S/N of 
the companion signal should be optimized. This optmization is 
typically achieved neither when the companion signal is not af¬ 
fected by any self-subtraction, nor when the speckle noise is best 
subtracted, leading to strong self-subtraction, but for some inter¬ 
mediate self-subtraction. 

To evaluate the impact of the algorithm throughput, we in¬ 
jected in the scientific datacubes simulated planets at different 
separations (0.2, 0.3, 0.4, 0.5, 0.6, 0.7 and 0.8 arcsec). To be 
able to clearly determine the optimal algorithm throughput for 
the single planets a wide J-band contrast of 3x10 was used. 
The contrast value was evaluated using the off-axis PSE not af¬ 
fected by the coronagraph. The spectrum of the simulated plan- 


Table 1. Summary of data used. 



YJ 

YH 

Date 

02/05/2013 

10/24/2013 

dit;^ 

4 s 

2s 

NDlf* 

56 

50 

Neutral Density Filter 

0.0 

1.0 

N datacubes 

16 

10 


“ Detector Integration Time 
* Number of frames for each datacube 


ets was assumed to be gray. In this case, we assume that the 
simulated companions do not have identical spectra to that of 
the central star. It is important to point out that the algorithm 
used to subtract the speckles makes use of spectral information, 
and the throughput is dep endent on the spec tral properties of the 
simulated planets (see e.g. lMaire et al.l20T4l for the case of dual 
band image). Moreover, the planet flux, which is quite bright in 
our case, could impact the spectral fitting. 

To simulate the rotation of the EOV, we inserted the sim¬ 
ulated planets with rotation of 5° between two successive dat¬ 
acubes leading to a full EOV simulated rotation of 75° for 16 
datacubes. The rotation considered here is of the order of 1° per 
minute corresponding at Paranal to the field rotation for a target 
with a declination between -10° and -38° and observed dur¬ 
ing its transit at the meridian. A continous rotation of 5° would 
of course introduce a smearing of a point source especially at 
the edge of the EOV. However, for this exercise we just wanted 
to check the speckle subtraction obtained from images taken at 
a distance of 5 minutes and a rotation of 5°. We, then, did not 
consider individual single exposures, rather we took a stack of 
the images at different times. On the other hand, in the real case 
the ADI algorithm should be modified to take the large smearing 
effect introduced by the rotation of the EOV into account. 

4. Results for IFS YJ-mode 

In this section we present the lES data reduction for its YJ-mode, 
using the procedures detailed in Section |3] An overview of the 
data used for this analysis is given in Table [T] 

4.1. Results obtained with the updated spectral 
deconvolution procedure 

We obtained very good speckle noise subtraction using the up¬ 
dated version of the SD and evaluating the attenuation of the 
satellite spots described above. We obtained attenuation factors 
ranging from 30 to 50, depending on the function used to repro¬ 
duce the speckle spectra. Thi s fac tor is similar to the results ob¬ 
tained by iTTiatte et al.l (l2007l) and ICrepp et al.l (l201 ill . However, 
to estimate the gain in limiting contrast, the algorithm through¬ 
put on the companion’s after the SD procedure must be con¬ 
sidered as described in Sectior lJAl The final reduced image ob¬ 
tained from the considered data set where simulated planets were 
inserted is displayed in figure|4] In this Eigure, the segmentation 
near the star center is generated by the normalization performed 
in different sectors as described in Section [34] In the external 
part of the image, the data reduction method is very efficient 
at removing these artifacts unlike at close separation where the 
method is less effective. 

In Eigure |5jwe show the 5cr contrast curves obtained by ap¬ 
plying the refined SD procedure described in Section I3.1I and 
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-144 -71 2 76 149 222 295 368 441 

Counts 

Fig. 4. Final reduced image obtained from the data set using the 
updated SD procedure where seven fake planets with contrast of 
3x10“^ were inserted. 


taking the algorithm throughput into account. In the following, 
the contrast curves are calculated using the median image from 
the final datacube produced by the speckle subtraction proce¬ 
dure. The exact pro cedure used to calcula te the contrast follows 
these steps (see e.g. lMasciadri et an2005h : 

- The standard deviation inside a 1.5AjD wide box centered 
on each pixel is calculated for each pixel of the image. 

- A median of these standard deviations is calculated for all the 
pixels located at the same separation from the central star. 

- The contrast at different separations is then calculated by di¬ 
viding these median values by the normalization factor. The 
normalization factor is obtained by multiplying the maxi¬ 
mum of the off-axis PSF, derived from the star not affected 
by the coronagraph, the ratio between the coronagraphic im¬ 
age single exposure and the off-axis image exposure, and the 
ratio of the transmission of the neutral density hlters used in 
the two cases to avoid the saturation of the detector. 

The same method was used in lEsnosito et al.l (1201 3h for the 
LET observations of HR8799 system giving good results when 
comparing the S/N in the hnal image with the S/N obtained from 
the contrast plot. As an example HR8799e was just above the 
detection limit in our contrast plot confirming the fact that it 
was clearly visible in the hnal image. Although it provides a 
good view of the contrast performance, this is not the most used 
method to calculate the contrast plot, which is usually obtained 
by deriving the standard deviation along a ring around the center 
of the image for each distance to the center. For this reason, we 
have dedicated further discussion to comparing the two methods 
in Appendix lAl 

Figure |5] shows the hnal results of this procedure. The con¬ 
trast solid curve in black is obtained from the original corona- 
graphic images, the green curve is obtained from the application 
of the SD and the orange and red dashed curves represent the 
contrast after the application of the simulated ADI procedure. 
These contrasts were obtained by applying the SD and htting 
a 1 degree function to obtain a good speckle subtraction while 
keeping the algorithm throughput to an high level. In this hgure 
we do not show the contrast curves corresponding to a separa¬ 
tion lower th an 0.2 arcsec (sma ller than the bifurcation radius 
as dehned bv lThatte et al.ll2()(y^ from the central star because at 



Separation (orcsec) 


Fig. 5.5 -cr contrast plot with the IFS operating in YJ mode while 
simulating observations of a bright star and using the updated SD 
method. The contrast solid curve in black is obtained without any 
speckle subtraction while the green curve is obtained applying 
the updated SD procedure. The red and orange dashed curves 
show the contrast obtained with the simulated ADI procedure in 
addition to the speckle subtraction. Different symbols are used 
to indicate the positions of known planets on this plot. 


these separations the method is no more effective. On this plot 
we added the positions of some known planets to enlighten the 
capability of IFS to retrieve their signal. Since the J magnitude 
for HD95086b is not known we have exploited its known H-K 
color-color diagrams to infer their possible magnitudes assum¬ 
ing the same J-Fl index. The error bar displayed in the figure 
account for the uncertainties of the method. We can see that we 
will be able to image ySPic b just using the SD and without apply¬ 
ing the ADI while the planets of HR8799 are close to the limit 
for speckle deconvolution alone. On the other hand, we will have 
to use the ADI to be able to image HD95086b. 

The algorithm throughput obtained for each planet is given 
in Figure |6] for htting function (used in the SD procedure) de¬ 
grees of 0 (corresponding to a constant value given by the me¬ 
dian of the spectra), 1, 2 or 3. In this hgure and in the following, 
the algorithm throughput has a value of 1 for no Hux loss of the 
fake planet while it has a value of 0 when the planet is com¬ 
pletely canceled out. Clearly the throughput is much larger with 
lower htting degrees. Besides, the difference between higher de¬ 
grees (2 and 3) is smaller than the difference between lower de¬ 
grees. Furthermore, as expected, the throughput is larger at larger 
separations. This occurs because companion objects extend on 
smaller fr actions of the speckl es held spectrum at large separa¬ 
tions (see lSnarks & For32002h . 

However, high-degree htting functions provide a better re¬ 
duction of the noise. For detection purposes, the best way to 
perform the SD is to optimize the S/N. These plots are shown 
in Figure |7] The S/N is obtained by dividing the planet signal 
integrated within a 1 pixel radius aperture centered on the planet 
position, by the standard deviation estimated in a 10x10 pixels 
region located near the considered planet but at a separation such 
that the Hux in this region is not affected by the planetary signal. 
This is an acceptable way to estimate the noise for our images 
since that the noise is very stable in the whole image. Figure Q 
shows that at separation lower than 0.35 arcsec and larger than 
0.75 arcsec, a ht with a low-degree polynomial function works 
better, while for intermediate separations the ht with a 3 degrees 
htting function gives the best S/N. Edge effects are caused by 
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Separation (orcsec) 


Separation (orcsec) 


Fig. 6. Algorithm throughput obtained for various degrees of the 
fitting function used in the SD procedure. A detailed explanation 
of the meaning of these curves is given in the text. 



Seporation (arcsec) 

Fig. 7. S/N plots obtained with different degrees of the htting 
function used in the SD procedure. A detailed explanation of the 
meaning of the plotted lines is given in the text. 


the noisy structure at the edge of the IFS FOV that can be seen, 
e.g., in Figure [3]or in Figure |4] and are important at separations 
greater than 0.6 arcsec. They can explain the steep fall of the S/N 
value at large separations obtained with high-order fits. 

Figure 0 displays the 5cr contrast curves obtained by taking 
the SD algorithm throughput (blue line) into account compared 
to the contrast obtained without taking it (red line) into account 
and before the application of the SD (green line). The difference 
from Figure |5] is that these contrasts are obtained just consider¬ 
ing the simulated planets. Despite this difference, the two latter 
curves are in good agreement with those shown in Figure |5] at 
least at separations larger than 0.3 arcsec. The discrepancy at 
0.3 arcsec can be given by the proximity to the central star and 
to the lower throughput of the SD at that separation. We have 
then tried to implement a very simple simultaneous spectral dif¬ 
ferential imaging (S-SDI) procedure by subtracting images ob¬ 
tained at different wavelengths after simulating a T7 spectral 
type planet. We have used the images at wavelengths between 
1.10 and 1.19 micron where the signal is not present to be sub¬ 
tracted from the images at wavelengths in the ranges between 
1.02 and 1.09 micron and between 1.22 and 1.29 micron where 
the signal from the planet is instead present. In this figure the 


Fig. 8. Contrast plots for the SD taking and not taking into ac¬ 
count the algorithm throughput. 


contrast curve obtained by applying this method is represented 
by the orange curve. The S-SDI is more effective at small separa¬ 
tions from the star (less than 0.3 arcsec) where the SD cannot be 
properly applied because the planet spectrum covers the entire 
spectru m provided by the IFS (see the definition of bifurcation 
point in iThatte et alfcOOTl) . 


4.1.1. Synthetic planets spectral retrieval 

We performed a test to check the ability of the system to re¬ 
trieve the spectra of the injected planets. The test was performed 
using the same data described in the previous paragraphs. The 
data were reduced using both the updated SD procedure alone 
and together with the simulated ADI procedure described in 
Section [33] The injected simulated planets had a mean contrast 
on the whole wavelength range of the IFS YJ-mode of 3x10“^. 
While the flux of the simulated planets was the same for all the 
instrument channels, the flux for the simulated star was variable 
and for this reason the contrast of the simulated planets changed 
for different channels. The output spectra for the planets at dif¬ 
ferent separations are shown in Figure|9]for the case in which we 
only used the SD procedure. In Figure[T0] however, we only used 
the ADI procedure. The latter gives a better reconstruction of the 
spectra and a lower dispersion of the results. In this second case 
all the spectra are very similar to each other indipendent from 
the separation from the central star. This is probably due to the 
high contrast of the simulated planets that we have used in this 
case. Flowever, in both cases we found a loss of flux with respect 
to the original spectra (represented by the black solid line in both 
the Figures). The final flux is around 44% of the original in the 
case of the spectral deconvolution alone, while it is of the order 
of 57% when we use the ADI procedure. This is due to the low 
algorithm throughput given by the speckle subtraction methods 
we describe. 

A more refined analysis of the spectro-photometric capa¬ 
bilities of the IFS and IRDIS and a description of the meth- 
ods adopted to im prove the algorithm throughput is reported in 
IZurlo et ^ (120141) . In that work the effects of different types of 
input spectra for the simulated planets are enlightened too. 
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Wavelength (micron) 


Fig. 9. Output spectra for the 7 simulated planets obtained ap¬ 
plying just the spectral deconvolution procedure. The black solid 
line represents the contrast for each IFS channel obtained from 
the input spectrum of the simulated objects. 



Wavelength (micron) 


Fig. 10. Output spectra for the 7 simulated planets obtained ap¬ 
plying the spectral deconvolution procedure and the simulated 
ADI procedure. The black solid line represents the contrast for 
each IFS channel obtained from the input spectrum of the simu¬ 
lated objects. 


4.2. Results with the PCA procedure 

Figure [TT] shows the hnal image obtained by applying the PCA 
procedure on a single science datacube without considering rota¬ 
tion of the FOV. Similar to the SD case, seven fake planets with 
3x10“^ contrast were injected at the same separation as above. 
In this hgure negative features along the direction star-planets 
arising from the subtraction procedure are clearly visible. They 
are caused by the oversubtraction given by a not perfect spectral 
htting around the planetary position. Figure[T2]displays the hnal 
image obtained using all the datacubes and where the rotation of 
the FOV is simulated with a rotation of 5° between consecutive 
datacubes corresponding to a total FOV rotation of 75°. Here we 
have decided to insert simulated planets with a lower contrast of 
10“^, which is nearer to the contrast limit for this image and has 
a smaller impact on a correct determination of the S/N probably 
because of the lower level of the strong subtraction structures 
described below. In the latter case, the image has been recon¬ 
structed using 32 principal components (see Section [3^ . while 
in the former case, without FOV rotation, only four principal 


0.1 arcsec 


•187 -112 -37 37 113 187 261 336 411 

Counts 

Fig. 11. Final image obtained by applying the PCA method to a 
dataset containing simulated planets with a contrast of 3x10“^. 
The image is obtained retaining only the hrst 4 principal compo¬ 
nents out of 38 and without considering the rotation of the FOV. 


• • • • • • • 


0.1 arcsec 


8.2 14 19 25 30 36 41 47 53 

Counts 

Fig. 12. Same as for Figure fTTIbut in this case the image is ob¬ 
tained retaining the hrst 32 principal components out of 608 and 
considering the rotation of the FOV too. Moreover, in this case 
the simulated planets have a contrast of 10^^. 


components were used. Negative structure given by the subtrac¬ 
tion algorithm are visible in Figure [121 The strong level of the 
noise suppression that we can obtain with this technique com¬ 
bined with the very high contrast of the simulated planets en¬ 
ables these structure to be visible. In future improvements of this 
procedure, we want to address in detail these residual structures 
in the image. However, the presence of these structures does not 
affect the scientihc results presented in this paper. Similar to the 
SD case, we exploited these images to calculate the algorithm 
throughput at different separations from the central star. The al¬ 
gorithm throughput versus the separation from the central star is 
shown in Figure [13] 

Using these algorithm throughputs, the contrast plot dis¬ 
played in Figure [14] were properly estimated. In this Figure the 
solid black line represents the contrast without any speckle sub¬ 
traction and the solid green line represents the contrast without 
held rotation (4 principal components out of 38 are used). The 
image resulting from the complete reduction of this data set is 
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Separation (orcsec) 


Fig. 13. Algorithm throughput versus separation for two sets of 
the PC A method parameters. The algorithm throughput is repre¬ 
sented by the orange curve and corresponds to no rotation of the 
FOV and 4 principal components, as shown in Figure [TT] while 
the green line represents the FOV rotation case and 32 principal 
components as shown in Figure [121 



Separation (orcsec) 


Fig. 14. 5cr contrast plot with the IFS operating in YJ mode for 
bright star observations after PCA data reduction and with simu¬ 
lated ADI procedure. The solid black line represents the contrast 
before speckle subtraction, the green line represents the contrast 
after PCA data reduction method using only the 4 principal com¬ 
ponents and without any rotation (see Figure [TTIl. The red line is 
obtained applying the PCA method with 32 principal compo¬ 
nents and with simulated FOV rotation (see Figure [T^. 


shown in Figure [TT] The dashed red line represents the contrast 
for the case where the FOV rotation is simulated and 32 princi¬ 
pal components out of 608 are used. The final image obtained 
from this data set is shown in Figure [12] 

We also focused on the variation of the S/N and of the algo¬ 
rithm throughput with the number of principal components used 
to reconstruct the original data. The throughput and S/N values 
for the planet at 0.4 arcsec separation are given in FigurefTSland 
in Figure [16] respectively. As expected, the throughput reduces 
with the number of principal components used, but at the same 
time the noise decreases so that the S/N is maximum when 32 
principal components are used. In order to find faint compan¬ 
ions, the number of principal components to be used needs to 
be carefully considered since it changes from case to case and. 


0 50 100 150 

Number of principal components 

Fig. 15. Algorithm throughput by the PCA method for a planet 
at 0.4 arcsec using different numbers of principal components. 


100.0 


10.0 

■Z. 

if) 

1.0 


0.1 

0 50 100 150 

Number of principol components 

Fig. 16. S/N for the PCA method with a planet at 0.4 arcsec using 
different numbers of principal components. 



for instance, it depe nds on the separation from the star (see e.g. 
iMeshkat et al .112014]) . However, for this analysis with simulated 
planets of known contrast we preferred to use a single number of 
principal components for the whole image. This choice is ade¬ 
quate for the objective of this analysis, which is mainly aimed at 
providing an initial estimate of the contrast that can be reached 
with this procedure. Moreover, this choice allowed us to use a 
simpler code and, as a consequence, to reduce the required com¬ 
puting time. 

Tests made using different numbers of principal components 
have allowed to us to find that it is adequate to use 32 of them 
as shown in Figure [14] We would like to stress that this num¬ 
ber is adequate just for the present data set. Different data sets 
will certainly require the use of a different number of principal 
components that must be evaluated on a case-by-case basis. 


5. Results for IFS YH-mode 

Following the scheme of Section ]4] we present the results of 
the analysis of the data taken with IFS operating in the YH- 
mode. As seen in Table [T] the data used for this instrument 
mode were taken at a different epoch with respect to the data 
taken for the YJ- mode. Between the two epochs the condition 
of the instrument changed. In particular the telescope simulator 
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Separation (orcsec) 


Fig. 17. Contrast plot obtained applying the SD and the ADI to 
the data taken with IFS operating in the YH-mode. The curves 
color codes are the same as in Figure |5] The difference in time 
delays given in the label with respect to those in Figure|5]is due 
to the different way in which the data have been taken at two 
epochs. See text for more detail. Different symbols are used to 
indicate positions of known planets on this plot. 


lamp was changed and an additional neutral density filter was 
used to avoid saturation. Moreover, the spectra obtained from 
the telescope simulator presented a strong absorbing band at 1.4 
yum caused by the OH ions contained in the optical fiber used. 
This caused extremely low flux at the wavelengths around 1.4 
yum making it more difficult to apply the SD. For these reasons 
the contrast obtained in this mode is worse than those for the 
YJ-mode. However, this problem should not occur on sky since 
the stellar spectra are flatter in this wavelength range with the 
resolving power of IFS, even considering the absorption due to 
water vapour in the Earth atmosphere. 

5.1. Results with the updated spectral deconvolution 
procedure 

In Figure [T7] the contrast plot obtained from the YH-mode 
IFS data is shown. These contrast levels are worse than those 
obtained with the instrument operating in the YJ-mode (see 
Figure|5]for comparison). The worsening is particularly evident 
for the simulated ADI contrast plot while it is much less evident 
for the SD alone contrast plot. This worsening could be due, be¬ 
sides the reasons given above, to the shorter total integration time 
for the data set used in this test. As for YJ-mode case, we have 
inserted the positions of some known planets on the plot. In this 
case, we will be able to image the planet of ySPicb without the ap¬ 
plication of the ADI, while to retrieve the signal from HD95086b 
and from the HR8799 planets, the use of ADI is mandatory. 

5.1.1. Synthetic planets spectral retrieval 

Similar to the YJ-mode analysis, seven simulated planets were 
added to the science datacube to estimate the precision in recon¬ 
structing the spectra of the injected planets. The introduced fake 
planets have the same flux for all the IFS channels. As a con¬ 
sequence, given that the flux of the simulated star was variable 
from one channel to the other, the contrast of the simulated plan¬ 
ets with respect to the simulated star was variable. The mean 
contrast on the whole wavelength range of the IFS YH-mode 



Wavelength (micron) 

Fig. 18. Spectral retrieval obtained for seven planets at different 
separations with the IFS operating in the YH-mode. In this case 
both the SD and the ADI were applied to the datacube before 
performing the spectral retrieval. The black solid line represents 
the contrast for each IFS channel obtained from the input spec¬ 
trum of the simulated objects. 



Wavelength (micron) 

Fig. 19. Spectral retrieval obtained for seven planets at differ¬ 
ent separations with the IFS operating in the YH-mode. In this 
case only ADI was applied to the datacube before performing 
the spectral retrieval. The black solid line represents the contrast 
for each IFS channel obtained from the input spectrum of the 
simulated objects. 


is of 3x10 In Figure [TSl the results obtained with this dat¬ 
acube after the application of both SD and ADI are shown, while 
Figure [19] shows the spectral retrieval obtained in this case after 
the application of ADI only on the original datacube. 

The spectral retrieval is better in the second case, confirming 
the results already obtained for the YJ-mode. At the same time, 
this result confirms that the application of SD strongly distorts 
the original spectra. This effect is slightly decreased when the 
ADI procedure is directly applied to the original datacube with¬ 
out applying the SD. Moreover, a second problem that affects 
the analysis for the YH mode is given by the great difference in 
flux between different images at different wavelengths in the YH 
datacube. This has the effect of inserting biases in the estimation 
of the normalization factors along the spectra and, in this way, 
leads to overestimate parts of the spectrum with respect to oth¬ 
ers. From the analysis of the results, we can conclude that the 
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Fig. 20. Contrast plot obtained applying the PCA procedure to 
the IFS YH-data. The solid green line represents the contrast 
obtained with the PCA applied to a single datacube without any 
rotation. The dashed red line represents the contrast obtained ap¬ 
plying the PCA to all the datacubes with the rotation of the FOV. 


Table 2. Summary of the contrast obtained at different separa¬ 
tions for the IFS YJ-mode. The meaning of the acronyms are 
explained in the text. 


Separation 

SD 

SD-tADI 

PCAl 

PCA2 

0.2 

1.33x10-'’ 

5.67x10-' 

1.30x10-'' 

2.51x10-® 

0.3 

5.94x10-® 

3.14x10-’ 

5.74x10-® 

1.11x10-® 

0.4 

4.60x10-® 

2.20x10-’ 

2.68x10-® 

6.59x10-’ 

0.5 

4.95x10-® 

2.36x10-’ 

2.51x10-® 

7.52x10-’ 

0.6 

6.31x10-® 

2.71x10-’ 

2.19x10-® 

8.20x10-’ 

0.7 

6.88x10-® 

2.90x10-’ 

1.45x10-® 

5.64x10-’ 

0.8 

1.93x10-® 

2.70x10-’ 

2.92x10-® 

4.05x10-’ 


distortion of the spectra is about eight times larger than in the 
case of the YJ-mode. Moreover, in some cases, additional flux 
is added to the original in such a way that the flux is higher in 
the extracted spectra than in the simulated spectra. This could 
arise from the worse quality of the data in the YH-mode than the 
YJ-mode as explained in Section lSTI Hence, careful analysis of 
YH-mode images will be then mandatory on on-sky data. 

5.2. Results with the PCA procedure 

Figure shows the contrast plots obtained applying the PCA 
to the IFS YH-mode data. The green line has been obtained ap¬ 
plying the procedure to a single datacube without any rotation 
of the FOV. We then applied the procedure to all ten datacubes 
to simulate FOV rotation as described in Section 14.11 The con¬ 
trast plot obtained for this case is represented by the dashed red 
line in Figure |20] As for all previous cases, the introduction of 
the field rotation allows us to obtain a contrast level better than 
lO-'’. 

6. Conclusions 

We have described the results of the functional and performance 
tests for the IFS instrument within SPHERE. While the func¬ 
tional tests performed using the instrument templates and the 
DRH recipes allowed us to validate and improve the software, 
the performance tests allowed us to get hints about the achiev- 


Table 3. Summary of the contrast obtained at different separa¬ 
tions for the lES YH-mode. The meaning of the acronyms are 
explained in the text. 


Separation 

SD 

SD+ADI 

PCAl 

PCA2 

0.2 

3.16x10-" 

3.32x10-® 

8.22x10-® 

1.90x10-® 

0.3 

1.02x10-" 

1.46x10-® 

6.04x10-® 

7.87x10-’ 

0.4 

6.86x10-® 

1.26x10-® 

1.01x10-® 

6.69x10-’ 

0.5 

6.60x10-® 

1.13x10-® 

9.95x10-® 

5.96x10-’ 

0.6 

6.00x10-® 

9.34x10-® 

7.45x10-® 

5.17x10-’ 

0.7 

4.73x10-® 

8.00x10-® 

5.47x10-® 

5.99x10-’ 

0.8 

3.11x10-® 

9.75x10-® 

1.43x10-® 

1.07x10-® 


able contrast with lES when used on sky and to optimize the 
calibration and data reduction procedures. The spectral decon¬ 
volution procedure currently implemented in the dedicated DRH 
recipe yields a contrast close to 10“^, w hile better results were 
expected from instrument simulations (iMesa et alj|201 Ih . Eor 
this reason we implemented an updated version of the SD pro¬ 
cedure taking several issues not properly addressed in the orig¬ 
inal procedure into account(e.g., more precise centering of the 
datacube images and more refined rescaling factor at different 
wavelengths). This updated procedure allowed us to improve the 
contrast by a factor of two to three. The final contrast obtained is 
of the order of 4x10“® at 0.3 arcsec separation, taking the total 
throughput of the differential imaging procedures properly into 
account. To consider this last effect, several simulated planets 
located at different separation from the star were injected. 

We also developed a new procedure for the speckle subtrac¬ 
tion exploiting the PCA algorithm. In Table |2] we summarize the 
5cr contrast we obtained with the lES in the YJ-mode with vari¬ 
ous speckle subtraction methods described in this paper. Similar 
results are presented in Table [3 for the lES YH-mode. In both 
these Tables the second column gives the results obtained when 
applying the SD to a single datacube, while in the third column 
lists the results obtained with our simulated ADI procedure af¬ 
ter SD are displayed. In the fourth column, labeled PCAl, we 
list the results obtained applying the PCA to a single datacube 
and with only four principal components used (in this case no 
rotation of the FOV was simulated), while in the fifth column, 
labeled PCA2, we give the results of the PCA performed on 
the complete data set in this case using 32 principal components 
(in this latter case a rotation of the FOV was simulated). While 
the cancellation of the noise is similar for both SD and PCA, 
the throughput is much lower for the PCA method as shown in 
Table |2l This latter effect seems to be much less important in 
the case of the YH-mode as shown in Table [3] As explained in 
Section|5] however, these data were taken under different condi¬ 
tions with respect to the YJ case and proper comparison between 
the two cases is difficult. Moreover this caused a loss of effective¬ 
ness of the SD method for the YH data (see Tabled. All these 
uncertainties will be solved when on-sky data become available. 

However, we can conclude that the PCA procedure can give 
results comparable to the SD results once the algorithm through¬ 
put is properly taken into account. A more refined method to 
take them into acc ount will be presented in an associated paper 
(IZurlo et al.ll2014ll . Moreover, the implementation of ADI pro¬ 
cedure in the PCA procedure is much easier than the implemen¬ 
tation of the ADI in the SD procedure. For this reason, the PCA 
procedure should probably be considered as a better solution for 
the final IFS pipeline. 

Finally, we have verified that, whatever the choice for the 
IFS reduction pipeline (SD, PCA or other possible solutions). 
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ADI is mandatory to reach the 10“^ 5cr contrast threshold using 
IFS mounted on SPHERE. The effectiveness of the ADI, how¬ 
ever, strongly depends on the rotation of the EOV and on the 
separation from the central star. 
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Appendix A: Contrast plot calculation 

The method that we have used to calculate the contrast plot in 
this paper and described in Section 14.11 has in the past been 
demonstrated to be reliable (see e.g. lEsposito et al.]|2013l) . On 
the other hand, the most common method used to calculate the 
contrast plot is to calculate the standard deviation on a ring cen¬ 
tered on the center of the star. Our choice was driven by the fact 
that this second method does not take the possible large scale 
inhomogeneities in the image properly into account, and that we 
need however to calculate some local noise when we are search¬ 
ing for a point source like a planet. 

To compare the results of the two methods we have run both 
of them on one of the cases presented here. In particular, we 
have chosen the case presented in Figure |5] On the left panel of 
Figure |AT| we present the same image while on the right panel 
we present the contrast obtained using the classical method de¬ 
scribed above on the same data set. The difference between the 
two plots (black line) for nonsubtracted images is readily evi¬ 
dent while it tends to diminish when we consider subtracted im¬ 
ages and it is almost null when we consider images where the 
ADI is applied. The great difference for nonsubtracted images is 
given by the problem of the large scale inhomogeneities in the 
original image cited above. To demonstrate this, we have pre¬ 
pared the same plots presented in Figure lATI but applying a low- 
frequency filtering on the nonsubtracted image, subtracting from 
the original image a smoothed version of the image itself using 
the SMOOTH IDL procedure with a large value of the smooth¬ 
ing area (with a semi amplitude of 25 pixels corresponding to 
~0.18 arcsec). 

The results of this test are presented in Figure IA.2I where 
it is evident that now the contrast plots are very similar using 
both methods. The conclusion of this test is that the main differ¬ 
ence between the two methods is on the nonsubtracted images 
where the large scale inhomogeneities in the image are not taken 
properly into account by the standard method. These differences 
can be reduced applying a low-frequency filtering on the images. 
The main problem with our method is that if the box that we 
are using for calculations has a dimension similar to the dimen¬ 
sion of the speckle (~ T/D), we could have an underestimation 
of the noise. However, we took care that the dimension of the 
box was larger than this value (as said in the paper, of the order 
of ~ 1.5/1/D). Moreover, various tests with different dimensions 
of the box have been done arriving at the conclusion that boxes 
with larger dimensions do not significantly change the estimated 
contrast. From these tests, we can estimate an error of the order 
of 10% for our contrast estimation. 

As a further test of the reliability of our contrast estima¬ 
tion we have introduced simulated planets with a contrast of 
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Fig.A.3. Nonsubtracted datacube median image with simulated 
planets with a contrast of 10“^. The planets with separation 
larger than 0.2 arcsec are just visible as foreseen by our contrast 
plot. 


10 on the nonsubtracted image (black plot on both panels of 
Figure lA.ll i. The result of this test is presented in Figure IA.3I 
where, as foreseen by the contrast plot, the planets at separations 
of 0.3 arcsec and larger are visible (even just above the detection 
limit) while the planet at 0.2 arcsec is not visible. This is in good 
agreement with what we can infer from the analysis of the con¬ 
trast plot displayed on the left panel and obtained exploiting our 
method. 

The conclusion of these tests are that both methods can give 
a reliable estimation of the contrast when used in the correct way. 
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Fig.A.l. Left: same of image |5] obtained with our method. Right: same contrast plot obtained using the standard method described 
in the text. 




Separation (orcsea) Seporotion (arasec) 

Fig.A.2. Left: same of image |5] obtained with our method where a smoothed version of the same image has been subtracted from 
the nonsubtracted image after applying a low-frequency filtering to subtract any possible large scale structure in the image. Right: 
same contrast plot obtained using the standard method described in the text. 























































